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MODEL FOR CALCULATING ELECTROLYTIC SHUNT PATH LOSSES IN 
LARGE ELECTROCHEMICAL ENERGY CONVERSION SYSTEMS 
by Paul R. Prokopius 
Lewis Research Center 

SUMMARY 

The quest for efficient energy conversion techniques may bring about the commer- 
cialization of electrochemical energy storage and generation systems. Within this gen- 
eral classification some of the systems may be designed with a common or circulating 
electrolyte supply. In such systems one of the factors that must be accounted for in the 
design is the power loss or self-discharge by cell to cell ionic shunt currents. This 
study addresses this design facet, the objective being to develop a mathematical analysis 
and solution technique to evaluate the shunt power losses in large electrochemical sys- 
tems. 

Anal 3 rtical tools consisting of an electrical shunt current circuit analog and a com- 
puterized solution technique were devised. This analytical work was done on the redox 
flow cell, but it is generally applicable in modified form to other systems. Representa- 
tive data are presented for a sample redox system, and the applicability of this analysis 
technique to other systems such as fuel cells or fluid flow networks is discussed. 


INTRODUCTION 

A design concept of many electrochemical energy conversion systems is to use a 
circulating or common electrolyte supply. Representative systems past and present in 
which this concept is employed are as follows: 

(1) Hydrazine - hydrogen peroxide fuel cells (Alsthom Company) 

(2) Low temperature alkaline fuel cell (Union Carbide) 

(3) Lithium - seawater battery (Lockheed Company) 

(4) Bulk energy storage redox (NASA- Lewis Research Center) 

The reasons for designing a system with a circulating electrolyte are varied. These 
include using the flowing electrolyte stream for a reactant delivery system as in the hy- 
drazine fuel cells and redox flow cells, a heat and water removal medium as in the 



alkaline cells, or a reactant control parameter as in the lithium -seawater battery. 
Though the individual design reasons vary, a parameter common to all systems which 
employ a circulating electrolyte is the power loss or self-discharge through cell to cell 
ionic shunt currents. 

This shunt loss characteristic needs to be considered in any system design; how- 
ever, most cases in the past did not require an indepth analysis. A common approach 
was to extrapolate basic shunt loss calculations to the size of the systems under consid- 
eration. In the simplest sense this extrapolation was made from sii^le shunt loop loss 
calculations or from small stack analyses. For the most part, this type of analysis was 
sufficient for the low voltage systems (i. e. , low shunt current driving voltage) built for 
imique applications - space, underwater, etc. In this type of system, while operating 
efficiency is an important consideration, design emphasis is placed on system reliability. 
This is not the case, however, for large terrestrial systems such as bulk storage redox. 
In this class of systems, design emphasis is placed on conversion efficiency as well as 
reliability; consequently, the shimt loss characteristic becomes more significant in the 
system design. 

With this background, a study was conducted to develop generalized shunt power loss 
analyses and solution techniques. An electrical circuit was proposed as an analog model, 
and a computer program was written with which to do the shunt loss calculations. The 
computer program was generalized in the sense that by making minor program input 
modifications virtually any size or type of electrochemical system can be analyzed. 


SHUNT CURRENT MODEL DEVELOPMENT 
Redox Concept 

A redox couple is characterized by a pair of oxidation-reduction reactions in which 
the ions of the couple remain soluble in their electrolytes in both the reduced and oxi- 
dized state. The rechargeable redox flow cell scheme is based on pumping two redox 
couple solutions through a power conversion unit to alternately charge and discharge the 
reactant solutions. In a utility bulk storage application the power conversion unit would 
be coupled with the power network such that the reactant solutions could be charged to 
store excess off-peak power and discharged to supply the peak power demand. A de- 
tailed description of the redox flow cell concept is presented in reference 1. 

A general schematic of a two-tank bulk energy storage redox system is shown in 
figure 1. The main components of this system are the catholyte and anolyte storage 
tanks, reactant circulating pumps, power conversion unit, and ac to dc and dc to ac 
power conditioning equipment. The couples chosen for illustration are titanium 
trichloride/ titanium tetrachloride and ferrous/ferric chloride. Their respective 
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charge-discharge reactions are shown on the schematic. In the discharge mode energy 
stored electrochemically in the electrol 5 de is converted and supplied on demand to the 
ac network to which the redox system is connected. Recharging is the reverse process 
and is accomplished by applying the appropriate voltage across the cell electrodes to 
drive current in the opposite direction. 

A redox system in a utility application would require a cell stack (power conversion 
unit) which operates in excess of 100 volts. A voltage level of this order of magnitude 
is needed to attain efficient operation of the high power semiconductors of the ac-dc 
power conditioners. The couples most commonly considered for application are aqueous 
solutions. Thus, the practical maximum attainable per cell charge voltage is about 
1. 2 volts because of the possible electrolysis of the water on charging. Therefore, in 
terms of stack size, a basic power unit of upwards of 100 series- connected cells would 
be required to produce the 100-volt minimum. Borrowing design characteristics de- 
veloped for fuel cells, the redox stack might, for the purposes of this study, be assumed 
to be assembled in a bipolar configuration; that is, the structural members separating 
the cells are conductive and provide the series electrical path through the stack. Also, 
a parallel reactant feed is assumed in which each of the individual cell reactant supply 
and return lines are paralleled off common manifolds. A schematic of a four-cell redox 
stack incorporating these two characteristics is shown in figure 2. 

Each individual cell of the power conversion unit consists of electrolyte supply 

chambers through which the reactant solutions are circulated, inert electrodes which 

serve as current collectors, and an anion- selective ion exchange membrane to separate 

the electrolytic solutions. In this configuration, the filled electrolyte supply/return 

lines connected by common manifolds form cell to cell shunt paths that are conductive 

to ionic currents. Ionic currents are generated and driven through these shunt paths by 

the cell to cell potential gradients of the stack. These currents represent an internal 
2 

I R power loss or self-discharge characteristic of the stack and thus must be kept to a 
minimum in any system design. The shunt losses become a system design parameter by 
virtue of the fact that the resistance of the shunt current paths is a function of the geom- 
etry of the electrolyte circulation system and the conductivity of the electrolyte 
(R„ = Z/GA where R_ is the shunt path resistance, I and A are shunt path length 
and area, and G is electrolytic conductivity). The tradeoff between shunt losses and 
electrolyte pumping power losses provide a simple example of one of the shunt loss de- 
sign considerations; that is, the shunt path resistance can be made very large and hence 
the shunt currents very small by designing the cell with long, small diameter electrolyte 
feed/return lines. On the other hand, an increase in shunt current resistance accom- 
plished in this manner is accompanied by an increase in electrolyte flow resistance 
which must be overcome by an increase in pumping power. To properly assess the shxmt 
loss portion of this design tradeoff as well as others involvii^ shunt currents, an ana- 
lytical technique is needed. 
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Circuit Analog 


For the system just proposed, a circuit analog of the stack current flow paths was 
devised as an analytical tool with which to do the shunt current analysis. In the analog 
each cell of the redox stack is represented by an ideal voltage source with an internal 
series resistance. The internal resistance is the electrical resistance through the cell 
in the direction of current flow. In each of the electrolyte shunt loops the electrical 
resistance to the ionic shunt currents is distributed along the path ler^th. However, in 
the analog circuit these resistances are lumped into passive resistor elements; that is, 
each of the four electrolyte supply and return taps of each cell is represented by a re- 
sistive element as is each manifold segment between any two taps. The complete analog 
circuit for the four-cell stack of figure 2 is shown in figure 3. The various components 
are represented as follows: cell internal resistance; e, ideal cell voltage (open 

circuit voltage); R^, anodic feed and exit port resistance; Rq, cathodic feed and exit 
port resistance; anode manifold segment resistance; cathode manifold seg- 

ment resistance; and Rj^, system load resistance. Several assumptions are inherent 
in the construction of the circuit analogy. First, it was assmned that the only compo- 
nent comprising the cell internal resistance R^ is the membrane resistance to ionic 
current flow. The resistance through the cell electrolyte in the direction of current 
flow (perpendicular to the electrode) plus that of the inert electrodes and the bipolar cell 
separator plates are small compared to the membrane resistance and are assumed 
negligible. Second, in keeping with the lumped parameter approach, the electrolyte in 
any of the reactant chambers is assumed to be at a uniform potential throughout; that 
is, the electrolyte discharge gradient from cell inlet to outlet is assumed negligible. As 
a result of this assumption, an internal phenomenon such as eddy currents generated by 
the discharge potential gradient across the cell is considered to be a second- order effect 
and is thus neglected. 

The analog circuit, by nature of the system being simulated, is a linear, symmetri- 
cal, passive, dc network to which Kirchoff's Law can be applied to determine the shunt 
current losses. Each cell except the two end cells of the stack generate four independ- 
ent conductive loops which have the cell equivalent circuit as a common element. The 
end cells generate only two loops by virtue of the fact that the conductive paths from the 
electrolyte supply to return manifolds are assumed broken in the pump circuit. Thus, 
the anolyte circuit of the cell on the low voltage end of the stack forms only one set of 
conductive loops with the anoljde circuit of the adjacent cell. Similarly, on the high 
voltage end of the stack, the cathol 5 d;e feed end return lines form only one set of loops 
with the adjacent cell. 
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MODEL ANALYSIS 


As stated previously, the shunt path analog circuit can be analyzed using Kirchoff's 
Law. The currents shown on the circuit diagram of figure 3 are the Kirchoff loop cur- 
rents for the 13 independent loops of the four-cell stack. The currents ij to ij ^2 
represent the shunt loop currents and ij^ the loop current for the load circuit, to this 
network, clockwise current flow was designated as the positive current direction. 
Equatii^ the sum of the IR voltage drops around each of the Kirchoff loops to the volt- 
age rise across the cell ideal source generates the voltage drop equations for the 13 in- 
dependent loops. These result in a set of 13 simultaneous algebraic equations to be 
solved for the 13 unknown loop currents. The loop currents thus identify the currents 
flowing through the various elements of the electrol 5 d;ic leakage paths, to some in- 
stances the loop current itself defines a leakage or shunt current factor. This occurs 
whenever the component in question is not in common with any other loop. The currents 
which are not determined directly are defined by a difference in loop currents such as 
(i 2 - i 0 ) which defines the leakage current flowing in the anolyte feed line to cell two. 

The leakage current through each individual resistive component is needed to calculate 
the power loss in each component. The individual component power losses can then be 
summed to determine the total shunt current power loss of the stack. 

For a redox stack of N cells with the stack configuration shown in figure 2, the ex- 
pression 4(N-1)+1 defines the number of Kirchoff loop equations in the circuit analogy. 
From this it is obvious that for any commercial sized stack, which would undoubtedly 
involve hxmdreds of equations, a computerized analysis is necessary. 

Of the many direct, iterative, or statistical methods of solving systems of linear 
equations. Grout's method appears to be well suited to this application in terms of com- 
putational and computer storage efficiency, to this method a matrix is constructed using 
the coefficients and independent variables of the Kirchoff equations. This coefficient 
matrix is manipulated on a term by term basis to obtain a derived matrix from which the 
equations' solutions can be obtained. The manipulation is such that as the various terms 
are being determined for the derived matrix, the terms of the original coefficient matrix 
need not be retained. Thus, only one matrix array need be allocated in the computer 
storage system. Also, because of the symmetry of the circuit being analyzed and the 
limited degree of commonality from loop to loop, the resultant equation matrix is largely 
symmetrical with only a few nonzero terms on either side of the major diagonal. Again, 
Grout's method lends itseK well in takii^ advantage of this symmetry and the large 
blocks of zero terms. The resultant effect is that a very large number of simultaneous 
equations can be handled with a very small array. 

The mechanics of the computational method may be best explained by an example. 
Given are three equations in the three unknowns, x^ X 2 , Xg: 
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^ 11^1 ^ 12*2 ^ ^ 13^3 ^ ^14 

^21^1 ^22^2 ^23^3 ^ ^24 

^31^1 '*' ^32^2 ^33^3 ^ ^34 

The coefficients and dependent variables of these equations form the matrix 

^11 ^12 ^13 ^14 

^21 ^22 ^23 ^24 

^31 ^32 ^33 ^34 

The various terms of this matrix are operated on to obtain a derived matrix 

^11 ^^12 ^13 -^14 

^21 ^22 ^23 ^24 

^31 ^32 ^33 -^34 

The expressions used in calculating the terms of the derived matrix are as follows: 

i-1 

^ii = ^ii " S ^iK^Ki 
K=1 

= i) (if i>j) 

K=1 

i-1 

Ay = - Z AikAk) (U 1 < i) 

K=1 J li 

Once the derived matrix is obtained, the solutions to the n unknown variables of the 
original equations are calculated from 
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n 


= E ^IK^K 

K=i+1 


[i = (n - 1)- ^ 


A computer program (Fortran IV) using Grout’s method was written for the redox 
configuration shov/n in figure 2. The program was written in a generalized form so that 
any sized stack of this particular configuration could be run directly or a variation from 
this configuration could be accommodated by a minimal modification in the program in- 
put relations. The actual stack size limitation is determined by the array storage limi- 
tation of the machine on which the program is being run. 

The program was written so that the only inputs required are the basic stack param- 
eters such as the number, of cells, the ideal cell voltage, and the various resistance 
values. From these data the equation and derived matrices are generated by the pro- 
gram and the loop currents are determined. In constructing the equation matrix, the 
two loops on either end of the stack plus the load circuit generate unique Kirchoff equa- 
tions; however, the interior loop equations along any one manifold are repetative, the 
only difference being an indexii^ of the current subscripts from loop to loop. Thus, the 
recursion equations for these interior loops were derived and programmed to facilitate 
generating the equation matrix. This matrix takes on a general form in which (1) all 
the terms within the body of the matrix are grouped in four locations to either side of the 
major diagonal, (2) all the positions in the last two columns are filled, and (3) the last 
row is completely filled. Read into a computer in its entirety, an equation matrix for 
just a 50-cell redox stack would occupy 38, 612 storage locations, a number which far 
exceeds the storage capacity of many machines. A characteristic of the Grout computa- 
tional technique is that the occurrence of zero in any location of the equation matrix will 
be duplicated by a zero term in the corresponding location of the derived matrix. This 
characteristic was used to expand the computational capacity of the program by elimi- 
nating all calculations involving zero terms of the equation matrix. The exception to 
this is the last row of the matrix in which each location is occupied by a nonzero term. 
This was handled as a special case in the logic by treating it as a single row matrix. In 
effect, the main body of the matrix was compressed from [4(no. of cells - 1) + 2] col- 
umns to only 11 columns in width. This width reduction logic was, by necessity, carried 
through the construction of the derived matrix and the computation of the results. 

To illustrate some of the general shunt current characteristics of a redox stack and 
to demonstrate the operation of the program, a hypothetical cell and stack design were 
assumed. The parameters that were used in this exercise are as follows: 

Electroljrte solutions: 

(Iron) 2M FeGlg in 0. 5N HGl 
(Titanium) 2M TiGlg in 6N HGl 
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Cell and stack geometry: 

2 

Active cell area, 2 ft 

Manifold to cell electrolyte feed and exit tubes, 1/8 in. i. d. , 6 in. long 
Supply and return manifolds, 3/4 in. i. d. 

Cell thickness, 0. 16 in. 

Membrane thickness, 0. 020 in. 

Stack size, 96 cells 

Supply and return manifold lengths, 15. 4 in. 

Electrical parameters: 

Iron couple specific conductance, 0. 151 (fi-cm)"^ 

Titanium couple specific conductance, 0. 165 (fi-cm)”^ 

Electrol 5 de feed and exit line resistance, 1205 
Cell to cell manifold resistance, 1. 0 
Membrane resistivity, 67 S2-cm 
Ideal cell voltage, 0. 7 V 
Load resistance, 0. 144 fi 

The system load resistance was set such that the stack would supply approximately 
200 amperes to the load circuit. For the anolyte feed manifold of a 96- cell stack, the 
shunt currents in the cell to cell manifold sections as calculated by the program varied 
from 7. 2 milliamperes for the end cells to 152 milliamperes in the center of the stack. 
The corresponding cell to cell voltage distribution varies from 0. 2890 volt for the end 
cells to 0. 2878 volt in the stack center. These distributions are shown in figure 4. The 
power delivered to the load by this 96-cell system is 6. 078 kilowatts, while the total 
shunt power loss was calculated to be only 11. 75 watts. This small shunt loss value was 
not unexpected since the electrolyte system was designed with a large feed line resist- 
ance. 

The system shunt power loss was found to increase exponentially with stack size. 
This characteristic is illustrated in figure 5 for the cell design assumed previously. The 
one component which is a tempering element in the slope of the exponential characteris- 
tic is the cell to cell manifold resistance since this represents an incremental resistance 
added to the circuit with each cell. 

The validity of the computer prc^ram was verified by comparirg results calculated 
with the program to those obtained from a standard simultaneous equation computer sub- 
routine. Because of the size limitation of the standard routine this comparison was done 
for a relatively small stack. However, because of the recursive natvire of the method 
programmed, validity can justifiably be assumed for any sized stack. 

CONCLUDING REMARKS 

The analysis technique described previously is not limited to the shunt current 
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analysis of the parallel feed, bipolar redox system. Examples of other possible uses 
include shimt current analyses of redox systems with other feed configurations or shunt 
current analyses of fuel cell systems employing circulating electrolytes. It could also 
be used for hydrodynamic analyses of flow networks such as the redox or fuel cell elec- 
trolyte supply systems. Each of these types of systems generate a modified form of 
the circuit analc^y that was developed for the two loop redox system. For example, the 
network analc^ for a fuel cell stack with a circulating electrolyte would consist of only 
two sets of Kirchoff loops instead of the four sets of the redox system. As another 

example, the network analog defining the fluid flow distribution of an electrolyte supply 

\ 

system would be represented by just a single set of Kirchoff loops. In this case the re- 
sistive elements of the redox analog would represent fluid flow resistances. The voltage 
sources would be omitted from the individual loops, but a voltage source representing 
the electrolyte pump pressure would be inserted in the load circuit. 

To accomplish a change in the programmed analysis from a four loop set to a two 
loop set network would merely require the omission of two sets of Kirchoff equations in 
constructii^ the equation matrix. Equally simple is the elimination or addition of ele- 
ments such as loop voltage sources since the equation matrix is constructed with the use 
of recursion equations. 

The general ease of program modification and the high degree of similarity between 
systems could be advantageous in doing system design studies. For example, the redox 
shunt current and electrolyte flow analyses could easily be combined in an iterative pro- 
gram for doing design optimization studies. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 15, 1975, 

506-23. 
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Rgure 1. - Two tank electrically rechargeable redox flow cell. 


Catholyte supply manifold 



Figure 2. - Four cell redox system flow schematic. 









